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SUMMARY 


An  analysis  of  the  steady-state  aerodynamic  heating  problem  at 
high-supersonic  speeds  is  made  for  two-dimensional  flows  with  laminar 
boundary  layers.  The  aerodynamic  heating  is  shown  to  be  reduced 
substantially  by  injecting  a  small  amount  of  coolant  through  a  porous 
surface  into  the  boundary  layer.  The  coolant  injection  acts  in  two 
ways  to  decrease  the  aerodynamic  heating:  First*  and  most  important, 
the  velocity  profile  is  altered  such  that  the  rate  at  which  heat  is 
conducted  to  the  surface  is  reduced  and,  second,  the  coolant  absorbs  an 
amount  of  heat  which  is  a  function  of  the  difference  in  temperature 
between  the  surface  and  the  coolant.  The  first  effect  provides  the 
advantage  of  cooling  by  injection  over  that  of  simply  using  a  coolant 
to  absorb  heat  from  the  surface.  Calculations  of  the  stability  of  the 
laminar  boundary  layer  show  that  for  a  wide  range  of  high-speed  flight 
conditions  the  boundary  layer  would  remain  laminar  at  all  Reynolds 
numbers  according  to  the  stability-theory  consideration.  The  analysis 
includes  calculations  of  the  cooling  requirements  and  equilibrium 
surface  temperatures  for  flat  plates  and  for  flat  porous  surfaces  with 
several  rates  of  fluid  injection  at  Mach  numbers  from  5  to  15  and  alti¬ 
tudes  from  sea  level  to  200,000  feet.  Some  calculations  of  the  skin 
friction  are  also  included. 


INTRODUCTION 


The  aerodynamic  heating  problem  assumes  considerable  importance  at 
high-supersonic  speeds.  Sanger  and  Bredt  (reference  l)  have  calculated 
the  high-speed  aerodynamic  forces  and  equilibrium  surface  temperatures 
at  extremely  high  altitudes  where  the  molecular  mean  free  path  is  large 
(free -molecule -flow  region)  compared  with  a  characteristic  body  dimension. 
Although  under  these  conditions  the  surface  temperatures  are  low,  the 
maximum  lift-drag  ratios  become  very  small.  Consequently,  flight  at  these 
altitudes  may  be  satisfactory  for  nonlifting  missiles  but  would  be 
unsatisfactory  for  steady  level  flight.  On  the  other  hand  in  the  lower 
atmosphere  (say,  sea  level  to  200,000  ft)  the  lift-drag  ratios  are 
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considerably  larger,  but  the  aerodynamic  heating  is  most  acute  and  some 
means  of  surface  cooling  must  be  employed  at  very  high  speeds.  The  scope 
of  the  present  paper  is  limited  to  a  consideration  of  the  steady-state 
aerodynamic  heating  problem  at  altitudes  where  the  air  may  be  considered 
a  continuum,  and  particular  emphasis  is  placed  on  a  means  of  cooling 
the  surface  by  continuous  injection  of  a  fluid  in  the  boundary  layer. 

Both  the  heat  transfer  and  the  drag  coefficients  are  known  to  be 
lower  for  laminar  than  for  turbulent  flows.  At  supersonic  speeds  it 
may  be  possible  to  maintain  a  laminar  boundary  layer  and  thus  alleviate 
the  aerodynamic  heating  somewhat.  The  theoretical  investigation  of  Lees 
(reference  2)  on  the  stability  of  the  laminar  boundary  layer  in  compres¬ 
sible  flow  indicates  that  the  laminar  layer  is  completely  stable  at  all 
Reynolds  numbers  at  supersonic  speeds  for  a  sufficiently  low  ratio  of 
surface  temperature  to  stream  temperature .  Thus,  for  flow  over  smooth 
surfaces  small  disturbances  are  damped  and  a  turbulent  boundary  layer 
does  not  develop  even  at  high  Reynolds  numbers.  The  possibility  of 
maintaining  a  laminar  flow  under  these  conditions  with  finite  but  small 
disturbances  in  the  boundary  layer  -  that  is,  maintaining  a  laminar 
boundary  layer  over  surfaces  that  may  be  employed  on  a  high-speed 
aircraft  -  has  not  been  verified  experimentally.  Experimental  verifica¬ 
tion  of  this  question  as  well  as  the  general  conclusions  of  Lees*  work 
is  therefore  desirable.  Nevertheless,  in  view  of  the  conclusion  of 
laminar  stability  for  infinitesimal  disturbances,  it  is  desirable,  to 
investigate  theoretically  methods  of  decreasing  the  aerodynamic  heating 
for  laminar  flows. 

The  aerodynamic  heating  is  alleviated  somewhat  with  increasing 
altitude,  but  even  with  laminar  flow  at  altitudes  of  200,000  feet  the 
surface  temperatures  may  be  excessive  at  high- supersonic  speeds.  An 
investigation  of  means  of  cooling  the  surface  is  therefore  desirable. 

One  possible  method  of  cooling  is  that  of  injection  of  a  cool  fluid 
through  a  porous  surface  into  the  boundary  layer  (reference  3)*  This 
coolant  injection  acts  in  two  ways  to  decrease  the  aerodynamic  heating: 
First,  it  alters  the  velocity  profile  such  that  the  heat-transfer  rate 
from  the  fluid  to  the  surface  is  reduced  and,  second,  the  coolant  absorbs 
an  amount  of  heat  which  is  a  function  of  the  difference  in  temperature 
between  the  surface  and  the  coolant.  The  first  effect  provides  the 
advantage  of  cooling  by  fluid  injection  over  that  of  simply  using  a 
coolant  to  absorb  heat  from  the  surface. 

Injection  of  the  coolant  through  a  porous  surface  affects  the 
boundary-layer  stability  in  two  ways.  The  direct  effect  of  injecting  a 
fluid  in  the  boundary  layer  is  to  alter  the  velocity  profile  such  that 
the  flow  is  less  stable.  The  indirect  effect  is  to  give  a  lower  surface 
temperature  which  in  turn  tends  to  make  the  flow  more  stable.  Actual 
calculations  must  determine  whether  the  flow  would  be  stable  for  given 
conditions.  This  question  will  be  discussed  in  more  detail.  The  blowing 
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rate  should  be  kept  as  low  as  possible,  consistent  with  adequate  cooling, 
since  the  coolant  will  probably  be  carried  by  the  aircraft.  Consideration 
of  storage  would  probably  dictate  the  use  of  a  liquid,  and  a  coolant 
which  has  a  low  temperature  and  a  high  beat  capacity  would  be  most 
efficient. 

In  the  present  paper,  calculations  of  the  steady-state  heating  of 
laminar  flows  are  presented  for  several  altitudes  from  sea  level 
to  200,000  feet)  Mach  numbers  of  5,  10,  and  15)  and  several  different 
surface  lengths.  Results  are  presented  for  a  flat  porous  surface  with 
three  rates  of  coolant  injection  as  well  as  for  a  flat  plate  (no  blowing) . 
For  the  case  of  continuous  injection  (injection  through  a  porous  surface) 
the  calculations  were  based  on  a  theory  for  Prandtl  number  of  1.0,  and 
in  all  cases  the  results  apply  to  a  flat  surface  at  zero  incidence. 
Included  in  the  present  paper  are  calculations  of  the  stability  of  the 
laminar  layer  and  some  calculations  of  the  skin  friction. 


SYMBOLS 


H 


c 


Hr 


Ha 

Hv 

AH 
T  . 
x 


average  heat  transfer  rate  at  surface,  Btu/(sq  ft) (sec) 

rate  at  which  heat  is  radiated  from  surface, 

Btu/(sq  ft) (sec) 

rate  at  which  heat  is  absorbed  by  coolant,  Btu/(sq  ft) (sec) 
heat  rate  to  vaporize  coolant,  Btu/(sq  ft) (sec) 

♦ 

heat-removal  rate,  Btu/(sq  ft) (sec) 
temperature,  °F  absolute 
length  of  surface  from  leading  edge,  feet 
coordinate  normal  to  surface,  feet 


k  conductivity,  Btu/(sec)(ft) (°F  abs. ) 

Cp  heat  capacity  at  constant  pressure,  Btu/(slug) (°F  abs.) 

h  coefficient  of  viscosity,  slugs/(ft) (sec) 
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P 

L 

a 


u 

v 


mass  density,  slugs  per  cubic  foot 
latent  heat  of  vaporization,  Btu/slug 
Stefan-Boltzmann  constant 

^4.8  X  10-13  Btu/(sec)(sq  ft)(°F  abs. 
velocity  parallel  to  surface,  feet  per  second 
velocity  normal  to  surface,  feet  per  second 


w 


u_ 

U1 


7 

M 

R 


Pr 


ratio  of  heat  capacities 
stream  Mach  number 


Reynolds  number 

Prandtl  number 


So 

5 


coolant  parameter 


2PgVg  ,  pjUjXX 

piui  V  h  J 


nondlmensional  stream  function 


D 

CDf 


skin  friction  drag,  pounds 
skin-friction  drag  coefficient 


Subscripts: 

s  value  at  surface 

i  coolant  value 

1  stream  value 

iso  isothermal  value 


# 


« 


* 
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t  stagnation  value 

re  recovery  value 


ANALYSIS 


In  order  to  determine  the  equilibrium  surface  temperature  of  a  body 
in  motion  through  a  fluid,  a  balance  between  the  heat  brought  to  and 
that  removed  from  the  surface  must  be  made.  This  transfer  of  energy 
occurs  mainly  through  the  processes  of  convection  and  radiation.  The 
analysis  is  simplified  by  considering  only  nighttime  operation.  Thus 
the  effects  of  solar  and  atmospheric  radiation  are  neglected  and  only 
radiation  from  the  aircraft  surface  is  considered.  With  these  restric¬ 
tions,  then,  the  heat  balance  may  be  written  as 


H  +  H  +  AH  =  0 
c  r 


(1) 


where  Hc  is  the  average  heat- transfer  rate  between  the  fluid  and  the 
surface,  Hr  is  the  rate  at  which  energy  is  radiated  from  the  surface, 
and  AH  is  the  rate  at  which  heat  is  absorbed  by  the  surface.  The 
temperature  at  which  AH  =  0  is  termed  the  natural  equilibrium  temperature. 
The  radiant  energy  for  a  grey  body  is  readily,  calculated  for  a  given 
surface  temperature  Ts  (°F  abs.)  from  the  expression 

Hr  =  0.9aTs^  (2) 


where  a  is  the  Stefan-Boltzmann  constant  and  the  factor  0.9  is 
chosen  as  the  emissivity  for  the  present  calculations.  The  rate  of  heat 
radiation  is  given  as  a  function  of  surface  temperature  in  figure  1.  The 
determination  of  the  conductive  heat  rate  requires  the  solution  of  the 
boundary -layer  equations.  The  numerical  solutions  for  high-speed  flows 
are  lengthy  and  tedious,  and  in  order  to  expedite  the  work  approximate 
solutions  were  used.  Although  this  procedure  affects  the  absolute  accuracy 
somewhat,  the  comparison  between  the  flat  plate  and  a  flat  surface  with 
fluid  injection  is  unaffected.  The  solutions  for  the  flat  plate  and  the 
surface  with  coolant  injection  differ  in  some  respects  and  are,  therefore, 
discussed  separately.  Specifically,  the  steady-state  heat  transfer  at 
supersonic  speeds  for  a  flat  plate  and  for  a  flat  porous  surface  with 
fluid  injection,  both  at  a  constant  temperature  over  the  entire  surface, 
are  considered  herein. 
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Conductive  heat  transfer.-  Throughout  the  analysis  the  flow  is 
considered  laminar,  and  the  temperatures  thus  calculated  are  later 
compared  with  the  results  of  the  stability  theory  of  the  laminar  boundary 
layer  to  determine  their  validity.  The  average  heat  rate  between  the 
fluid  and  a  surface  of  unit  width  and  of  length  x  may  be  written 
according  to  Fourier !s  law  of  heat  transfer  as 


where  k  is  the  conductivity  of  the  fluid,  T  is  the  temperature, 
y  is  the  coordinate  normal  to  the  surface,  and  the  subscript  s  denotes 
the  value  at  the  surface.  The  quantity  dT/dy  is  found  by  solving  the 
dynamic  and  energy  equations  of  the  boundary  layer.  Even  for  a  simple 
configuration  such  as  a  flat  plate  the  problem  as  stated  becomes  very 
difficult  for  high-speed  gas  flows.  Several  investigators  have  calculated 
velocity  and  temperature  profiles  from  the  boundary-layer  equations  by 
various  numerical  integration  methods  (reference  4,  for  example)}  however, 
these  calculations  are  insufficient  for  the  present  analysis.  The  usual 
method  of  determining  the  high-speed  heat  transfer  (references  5  and  6) 
is  to  modify  the  low-speed  solutions.  This  practice  is  followed  herein. 

Conductive  heat  transfer  -  flat  plate. -  The  usual  generalization 
of  the  low-speed  heat-transfer  equations  results  by  employing  values  of 
density,  viscosity,  conductivity,  and  heat  capacity  based  on  surface 
conditions  instead  of  stream  values  and  by  modifying  the  temperature 
potential  (reference  5)*  For  laminar  flow  over  a  flat  plate  then  the 
average  heat-transfer  rate  between  the  fluid  and  surface  is 


where  the  subscript  0  denotes  the  value  for  a  flat  plate,  pQ  is  the 

density  at  the  surface,  [±q  is  the  viscosity  at  the  surface,  Pr  is  the 
Prandtl  number  (taken  as  0.72),  U]_  is  the  stream  velocity,  and  Tre  is 
the  recovery  (or  adiabatic  wall)  temperature,  that  is,  the  temperature 
an  insulated  nonradiating  surface  would  assume  in  an  adiabatic  flow.  This 
form  of  the  heat-transfer  equation  has  proved  satisfactory  at  moderate- 
supersonic  speeds,  however,  it  has  not  been  verified  experimentally  at 
the  high  speeds  considered  herein. 
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The  ratio  of  the  temperature  rise  with  no  heat  transfer  to  the  stagnation- 
temperature  rise  is  termed  the  recovery  factor.  A  knowledge  of  this 
factor  permits  the  determination  of  the  recovery  temperature.  In  general 
this  factor  is  dependent  upon  the  geometry  of  the  body  and  the  Prandtl 
number  and  is  found  from  experiment  or  calculation.  A  value  of  0.9 
is  used  herein  and  the  recovery  temperature  then  is  found  from  the 
relation 


Tre  =  T-l  +  0.9 (Tt  -  Ti)  (5) 

where  is  the  stream  temperature  and.  is  the  stagnation  tempera¬ 

ture.  In  order  to  allow  for  the  variation  in  the  ratio  of  specific 
heats  with  temperature,  the  stagnation  temperature  was  calculated  from  a 
relation  derived  by  quantum  statistical  methods  (reference  6) 


where  M  is  the  stream  Mach  number.  Figures  2  and  3  present  the 
variation  of  stagnation  and  recovery  temperatures,  respectively,  with 
Mach  number  for  several  altitudes. 

The  heat  balance  for  laminar  flow  over  a  flat  plate  at  several 
altitudes  and  Mach  numbers  was  determined  from  equation  (l).  The  altitude 
tables  of  reference  7  were  used  to  determine  the  stream  temperature  and 
density  and  the  viscosity  and  conductivity  coefficients  were  calculated 
from  Sutherland's  equation.  The  heat  absorbed  by  the  surface  (AH)  for 
a  10-foot  surface  length  is  given  as  a  function  of  the  surface  tempera¬ 
ture  in  figure  k  for-  several  altitudes.  The  natural  equilibrium 
temperature  is  found  from  the  condition  that  AH  =  0.  This  temperature, 
for  the  flat  plate  (£0  =0),  is  given  in  figure  5  for  several  altitudes, 
Mach  numbers,  and  surface  lengths.  In  general  the  equilibrium  temperature 
is  found  graphically  or  by  trial  and  error. 

Heat  balance  with  continuous  fluid  injection.-  The  high  natural 
equilibrium  temperatures  of  a  flat  plate  at  high- super sonic  speeds 
emphasize  the  need  for  some  adequate  means  to  cool  the  surface.  Injection 
of  a  fluid  through  a  porous  surface  into  the  boundary  layer  changes  the 

chi 

velocity  profile  such  that  the  velocity  gradient  —  at  the  surface 
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is  decreased.  Since  the  temperature  gradient,  and  thus  the  heat-transfer 
rate,  is  approximately  proportional  to  the  velocity  gradient,  the  rate  at 
which  heat  is  conducted  to  the  surface  is  reduced.  A  very  small  rate  of 
fluid  injection  decreases  the  aerodynamic  heating  significantly.  Since 
the  quantity  of  fluid  injected  is  small,  the  usual  boundary-layer  equations 
apply. 

If  the  initial  temperature  of  the  fluid  injected  into  the  boundary 
layer  is  lower  than  the  surface  temperature,  the  coolant  will  absorb 
an  amount  of  heat  proportional  to  the  temperature  difference j  furthermore, 
if  the  coolant  is  a  liquid,  it  will  absorb  an  additional  amount  of  heat 
proportional  to  the  heat  of  vaporization.  '  The  change  in  the  velocity 
profile  also  alters  the  stability  of  the  boundary  layer.  In  reference  3 
it  is  shown  that  the  fluid  injection  is  destabilizing,  but  on  the  other 
hand  the  effect  of  the  heat  transfer  is  stabilizing.  Thus,  the  net. 
effect  of  injection  and  cooling  must  be  accounted  for  in  the  stability 
consideration.  For  a  wide  range  of  conditions  the  flow  is  stable  in  the 
laminar  form.  See  section  entitled  "Stability  considerations  of  the 
laminar  boundary  layer." 

For  the  present  analysis  it  is  assumed  that  the  coolant  is  liquid 
air  and  that  part  of  the  heat  in  the  boundary  layer  is  utilized  to 
vaporize  the  coolant.  The  heat  balance  then  may  be  written  as 


Hc  +  Ha  +  Hy  +  Hr  =  0 


(7) 


where  Ha  is  the  rate  at  which  heat  is  absorbed  by  the  coolant  and  Hv 
is  the  heat  rate  required  to  vaporize  the  coolant.  The  heat  radiated 
from  the  porous  surface  is  given  by  equation  (2)  and  the  average  value 
of  Hv  is 


x  ,  ?oplulL 
LPsvs  _  2x 


UiPfX 


dx-  =  -l£ 


Plul 

\T? 


(8) 


where  L  is  the  latent  heat  of  vaporization  (taken  as  283^  Btu/slug), 

R  is  the  Reynolds  number  based  on  stream  condition,  pg  is  the  density 

of  the  coolant,  v3  is  the  velocity  normal  to  the  sin-face,  and  £0  is 
a  nondimensional  flow  parameter  defined  as 


bo 


2Psys  ftP 1* 

piui  V 


(9) 
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The  parameter  £0  arises  in  the  treatment  of  the  dynamic  and  energy 

equations  (reference  3)«  The  condition  of  constant  temperature  over  the 
surface  is  satisfied  for  =  Constant. 

In  a  similar  manner  the  average  rate  at  which  heat  is  absorbed  by 
the  coolant  is  found  as 


Ha  =  x  psvscPi(Ts  -  Ti)to 


P 1  "Ll-i 

Ha  =  -  Ti) 


where  Cp^  is  the  specific  heat  and  is  the  temperature  of  the  vaporized 

coolant.  The  heat  capacity  cpi  was  taken  as  7-712  Btu/(slug)(°F  abs.) 

and  the  temperature  T-^,  as  1^7°  F  absolute. 

As  in  the  case  of  the  flat  plate,  the  calculation  of  the  heat 
transfer  between  the  boundary  layer  and  the  surface  requires  a  solution 
of  the  boundary-layer  equations.  However,  for  the  case  of  continuous 
injection  a  solution  with  certain  simplifying  assumptions  is  readily 
determined  from  known  isothermal  solutions  by  a  simple  quadrature 

(reference  3).  Then  the  quantity  (  — \  may  be  evaluated  in  terms  of 

Ws 

the  slope  of  the  velocity  profile  for  the  isothermal  problem,  and  the 
average  heat-transfer  rate  between  the  boundary  layer  and  surface  is 
readily  found  from  equation  (3). 

With  the  assumption  that  the  Prandtl  number  is  unity,  the  dynamic 
and  energy  equations  have  the  same  integral  in  the  absence  of  a  pressure 
gradient  and  the  temperature  is  related  to  the  velocity  by 


1  =  P-i  =  is  -  !s  .  fi  +  w 
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where  w  is  defined  as  u/u]_  and  7  is  the  ratio  of  .specific  heats. 

Reference  3  shows  that,  for  flow  over  a  porous  surface  of  uniform 
temperature  with  continuous  injection,  the  equations  of  the  compressible 
flow  are  of  the  same  form  as  the  isothermal  case  provided  the  Prandtl 
number  is  unity  and  the  viscosity  varies  linearly  with  temperature .  The 
coordinate  q  then  is  related  to  the  isothermal  value  by 


where  the  subscript  iso  refers  to  the  isothermal  value  and  C  is  a 
nondimensional  stream  function.  From  equations  (ll)  and  (12)  the 


quantity  /  — )  (the  value  at  w  =  0)  can  be  expressed  in  terms  of  the 

Ws 

isothermal  values  of  the  slope  of  the  velocity  profile  and  the  average 
convective  heat  rate  is  then  found  as 


(13a) 


or 


Hc  =  (Hc)5 


is  tabulated  for  several  values 

of  £  in  reference  8.  The  quantity  (Hc)  is  the  convective  heat 

o  b 

rate  for  no  blowing  -  that  is,  for  a  flat  plate. 


A  quantity  equal  to  four  times 


/  dw 

\dq 


dw\ 

®i/=isO 


„=o  r 


_ 

dl1/siso., 

-^o=0 


(13b) 


This  relation,  however, 
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does  not  give  exact  agreement  with  the  results  previously  given  for  the 
flat  plate.  In  order  to  provide  a  better  basis  for  comparison  of  the 
equilibrium  temperatures  with  and  without  blowing,  the  quantity  (Hc)«. 

'  *  o 


=0 


may  be  replaced  by  HCq;  therefore,  the  convective  heat  rate  with  blowing 
may  be  written  as 


The  heat  balance  for  surfaces  with  coolant  injection  is  determined 
from  equation  (7)  by  employing  the  foregoing  relations.  The  equilibrium 
temperatures  for  several  blowing  rates,  altitudes,  and  Mach  numbers  are 
given  as  a  function  of  surface  length  in  figure  5. 

Throughout  this  paper,  the  boundary  layer  is  assumed  to  be  laminar, 
and  this  condition  may  not  be  fulfilled  for  some  temperatures  calculated. 
The  results  should  then  be  checked  against  the  stability  criterions  for 
laminar  boundary  layers. 

Stability  considerations  of  the  laminar  boundary  layer.-  The  direct 
effect  of  injection  of  a  fluid  in  the  boundary  layer  is  to  destabilize 
the  flow*  whereas,  the  effect  of  removing  heat  is  stabilizing.  The 
stability  problem  then  is  to  determine  the  minimum  critical  Reynolds 
number  (that  is  the  Reynolds  number  at  which  the  flow  first  becomes 
unstable)  for  a  given  surface  temperature,  Mach  number,  and  amount  of 
coolant  injection.  The  range  in  which  the  heat-transfer  calculations  are 
valid  can  therefore  be  determined  by  comparing  them  with  the  results  of 
such  stability  calculations.  In  reference  2,  a  detailed  consideration  of 
the  stability  problem  is  presented  for  compressible  viscous  fluids,  and  in 
reference  3*  an  analysis  is  presented  for  the  specific  case  of  flows  over 
a  porous  surface  with  continuous  injection. 

Instead  of  determining  the  minimum  critical  Reynolds  number,  it  is 
more  convenient  to  determine  the  surface  temperature  required  to  make  the 
flow  stable  at  all  Reynolds  numbers  for  a  given  Mach  number  and  fluid 
injection  rate.  Because  of  the  extreme  curvature  of  the  velocity  profiles 
with  injection  of  gas  at  the  wall,  the  approximate  formulas  of  reference  2 
for  the  minimum  critical  Reynolds  number  are  not  applicable  and  the 
iteration  method  of  that  paper  must  be  employed.  An  estimate  of  the  value 
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of  minimum  critical  Reynolds  number,  which  will  serve  as  a  stability 
criterion,  is  obtained  by  taking  the  phase  velocity  to  have  its  maximum 
possible  value.  This  condition  occurs  when  ^(z)  =  O.580.  (See  refer¬ 
ence  2  for  the  definition  of  the  stability  parameters.)  For  that  value 
of  $i(z),  the  two  branches  of  the  curve  of  wave  number  a  against  R 
coincide  and  hence  the  phase  velocity  for  each  branch  approaches  the 

same  value  (c  =  1  -  — 

\  M 

The  actual  calculation  consists  of  choosing  a  stream  Mach  number 

and  fluid  injection  rate.  The  value  of  c  is  taken  as  1  -  —  and  a 

M 

trial  value  of  Ts/Tp  is  assumed.  Equations  (24)  and  (25)  of  refer¬ 
ence  2  are  used  to  obtain  a  zeroeth  approximation  to  the  values  of  $i(z) 

and  $r(z).  Figure  9  and  equations  (28)  and  (29)  of  reference  2  are  used 
to  determine  successive  approximations  to  ^ i ( z) .  If  the  final  value 
of  <E>i(z)  £  O.580  the  assumed  value  of  Ts/Tp  was  incorrect  and  a  new 
value  must  be  tried  with  the  corresponding  new  velocity  profile. 


as  R 


The  surface  temperature  for  stability  is  given  as  a  function  of 
Mach  number  for  a  given  amount  of  injection  in  figure  6,  and  those 
temperatures  may  be  compared  directly  with  those  found  from  the  heat- 
transfer  calculations.  If  the  calculated  surface  temperatures  are  equal 
to  or  less  than  those  found  from  the  stability  considerations  for  a 
given  Mach  number  and  coolant  injection  rate,  the  flow  is  stable  for  all 
Reynolds  numbers.  For  finite  values  of  the  minimum  critical  Reynolds 
number  somewhat  higher  surface  temperatures  for  stability  would  be 
allowed. 

Skin  friction.-  The  skin  friction  is  readily  determined  once  the 
velocity  profile  is  known.  The  skin  friction  drag  is  given  by 


■  * 
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or,  in  terms  of  drag  coefficient, 


(15) 


For  a  linear  variation  of  viscosity  with  temperature,  equation  (15) 
reduces  to 


D 


2"Plul 


4  Mb  ^s/dv^ 

\[R  Ml  Pl\dT)y 


Iso 


4  /dw\ 

>/rU/8. 

V  slSO 


(16) 


The  skin-friction  drag  coefficient  calculated  from  equation  (l 6)  is 
presented  in  figure  7*  This  quantity  is  independent  of  Mach  number 
because  of  the  linear  viscosity-temperature  relationship.  The  skin- 
friction  coefficient  decreases  with  Mach  number  (Prandtl  number  less 

than  1.0)  for  a  viscosity- temperature  relation  given  by  — 

Ml 

(See  reference  4.)  The  actual  variation  with  Mach  number,  however, 
depends  to  a  large  extent  on  the  nature  of  the  viscosity-temperature 
relation.  In  general  the  variation  of  the  drag  coefficient  with  Mach 
number  is  small  compared  to  the  change  due  to  fluid  injection. 


rT  \O.76 


DISCUSSION 


The  results  of  the  steady-state  heating  analysis  for  laminar 
boundary  layers  are  shown  in  figures  4  and  5.  Although  the  analysis 
refers  specifically  to  a  flat  surface  at  zero  incidence,  the  calculations 
are  indicative  of  the  high-speed  heating  problem  in  general.  The 
calculations  show  that  the  steady-state  temperatures  are  significantly 
dependent  on  the  Mach  number  and  altitude  and  to  a  lesser  extent  on  the 
surface  length,  and  that  continuous  injection  of  a  cool  fluid  in  the 
boundary  layer  reduces  the  equilibrium  temperatures  substantially. 

The  average  amount  of  heat  that  must  be  removed  from  the  surface 
of  a  flat  plate  10  feet  long  to  maintain  a  given  equilibrium  surface 
temperature  is  shown  in  figure  4.  It  is  readily  apparent  that  the 
cooling  requirements  are  increased  with  high  speeds  and  low  altitudes. 
Similar  calculations  for  other  surface  lengths  indicate  that  the  aero¬ 
dynamic  heating  is  alleviated  somewhat  for  longer  surfaces. 
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The  calculated  equilibrium  surface  temperatures  are  presented  in 
figure  5  as  a  function  of  the  length  of  the  surface  for  several  altitudes, 
Mach  numbers,  and  coolant  injection  rates.  It  is  apparent  that  the 
temperatures  of  the  surfaces  with  coolant  injection  are  appreciably 
lower  than  those  of  the  flat  plate  (Co  =  0) •  For  these  calculations  the 
coolant  was  assumed  to  be  liquid  air.  By  choosing  a  coolant  with  a  higher 
heat  capacity,  even  lower  temperatures  would  result.  From  consideration 
of  a  high  heat  capacity  and  a  low  vapor  temperature,  liquid  hydrogen  would 
be  a  logical  choice  for  the  coolant.  The  injection  of  foreign  gas  may 
alter  the  stability  considerations  somewhat,  and  for  that  reason  no 
calculations  for  the  case  in  which  hydrogen  is  used  as  the  coolant  are 
included.  Figures  4  and  5  indicate  that  the  equilibrium  temperatures 
increase  rapidly  with  Mach  number  and  are  highest  for  low  altitudes.  The 
effect  of  surface  length  on  the  equilibrium  temperatures  is  not.  very 
significant.  At  the  lower  values  of  the  coolant  parameter  Co  the  short 
lengths  are  somewhat  hotter,  but  at  the  larger  values  of  C0  there  is 
little  variation  with  length.  The  amount  of  coolant  injection  required 
to  provide  a  sizeable  decrease  in  the  convective  heat  rate  is  relatively 
small.  For  example,  for  the  most  extreme  condition  considered  here 
(lO -ft  surface  length,  sea  level,  M  =  15,  £0  =  _1),  t:iie  avera6e  mass 
flow  is  0.00122  slug/(sq  ft) (sec).  At  an  altitude  of  100,000  feet  this 
value  reduces  to  0.00012  slug/(sq  ft) (sec).  At  a  Mach  number  of  5  the 
mass  flow  is  0.0007  slug/(sq  ft) (sec)  and  0.00007  slug/(sq  ft) (sec)  at 
altitudes  of  sea  level  and  100,000  feet,  respectively.  For  a  flow 
parameter  of  Co  =  -0.5  the  mass  flow  would  be  one-half  of  those  given. 

The  average  coolant  mass  flow  decreases  for  increased  surface  lengths. 

The  primary  effect  of  coolant  injection  is  the  reduction  in  the 
amount  of  heat  transferred  from  the  boundary  layer  to  the  surface.  This 
effect  may  be  seen  from  equation  (l4).  The  heat-transfer  rate  with 
injection  is  equal  to  that  for  the  flat  plate  times  the  ratio  of  the 
slopes  of  the  velocity  profiles,  and  this  ratio  is  0.492,  0.282,  and  0.107 
for  Co  equal  to  -0.5,  -0.75,  and  -1.0,  respectively.  Thus,  the  heat 
rates  at  the  surface  for  Co  =  "0.5  and  -1.0  are  approximately  1/2  and  1/10, 
respectively,  of  that  of  the  flat  plate. 

The  fact  that  the  heat-transfer  rate  with  injection  is  appreciably 
lower  than  that  of  the  flat  plate  suggests  the  possibility  that  part  of 
the  boundary  layer  near  the  trailing  edge  could  be  sucked  off  and  used  as 
the  fluid  to  be  injected  in  the  boundary  layer.  This  air,  of  course,  is 
hot  and  would  heat  up  the  boundary  layer,  but  it  might  be  expected  that 
the  net  effect  would  be  a  reduction  in  surface  temperature  over  that  of 
a  flat  plate.  A  few  calculations  were  made  on  this  basis,  but  for  all 
conditions  calculated  the  heat  added  was  found  to  be  large  and  the  net 
effect  was  to  produce  equilibrium  surface  temperatures  slightly  higher 
than  those  for  a  flat  plate. 
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The  critical  ratio  of  surface  temperature  to  stream  temperature 
for  stability  of  the  boundary  layer  for  all  Reynolds  numbers  is  given  as 
a  function  of  the  stream  Mach  number  in  figure  6.  Accurate  calculation 
of  the  stability  at  high  Mach  numbers  is  difficult  and  for  that  reason 
the  results  are  shown  with  dashed  curves  above  a  Mach  number  of  5*  For 
all  the  results  in  figure  5,  however,  the  stability  analysis  indicates 
that  the  boundary  layer  is  laminar  for  all  Reynolds  numbers.  The 
effect  of  injection  on  the  stability  is  clearly  shown  in  figure  6.  At 
high-supersonic  speeds  (say  M  >  3)  the  surface  radiation  alone  provides 
adequate  cooling  to  maintain  a  laminar  boundary  layer  over  smooth 
surfaces . 

The  variation  of  skin-friction  drag  coefficient  with  Reynolds 
number  is  shown  in  figure  7*  The  fact  that  the  skin  friction  is 
independent  of  Mach  number  arises  from  the  assumption  of  a  linear 
viscosity-temperature  relationship.  For  other  viscosity-temperature 
relations  the  skin  friction  shows  some  variation  with  Mach  number,  and 
the  actual  variation  depends  on  the  nature  of  the  viscosity-temperature 
relation.  The  skin  friction  is  relatively  small  and  decreases  with 
increased  blowing  rates  as  long  as  the  boundary  layer  is  laminar. 


CONCLUDING  REMARKS 


An  analysis  of  the  steady-state  aerodynamic  heating  problem  at 
high-supersonic  speeds  for  two-dimensional  flows  with  laminar  boundary 
layers  shows  that  the  aerodynamic  heating  is  reduced  substantially  by 
injecting  a  small  amount  of  coolant  through  a  porous  surface  into  the 
boundary  layer.  Calculations  of  the  stability  of  the  laminar  layer 
show  that  for  moderate  rates  of  injection  the  boundary  layer  is  stable 
at  all  Reynolds  numbers,  and  these  injection  rates  provide  adequate 
cooling  for  many  high-speed  flight  conditions.  The  calculations  show 
that  the  equilibrium  surface  temperatures  are  significantly  dependent 
on  the  Mach  number  and  altitude  and  to  a  lesser  extent  on  the  surface 
length.  In  view  of  the  conclusion  of  laminar  stability  and  the  large 
decrease  in  aerodynamic  heating  with  a  small  amount  of  coolant 
injection,  it  is  desirable  to  investigate  the  problem  further  both 
experimentally  and  analytically  since  the  method  appears  practical 
for  some  high-speed  flight  applications.  Further  investigations  should 
include  a  study  of  the  effectiveness  of  injecting  a  fluid  into  the 
laminar  sublayer  to  decrease  turbulent  aerodynamic  heating. 
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Stagnation  temperature,  °F 
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Figure  2 Effect  of  Mach  number  on  stagnation  temperature. 


Recovery  temperature,  °P 
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Figure  3*“  Effect  of  Mach  number  on  recovery  temperature. 


Figure  4.-  Average  heat-removal  rate  for  flat  plate  of  10-foot  length 
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(d)  Altitude,  150,000  feet 


Figure  b.-  Concluded. 


Equilibrium  surface  temperature,  °F 
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Equilibrium  eurface  temperature, 
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Surface  length,  ft 

(b)  Altitude,  50,000  feet. 
Figure  5  *“  Continued. 


Equilibrium  surface  temperature,  °F 
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(c)  Altitude,  100,000  feet. 


Figure  5*-  Continued 


Equilibrium  surface  temperature, 
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Critical  surface  temperature  for  boundary  layers  that  are  stable  for  all  Reynolds  numbers 
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